Transport in finite size systems: an exit time approach 
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Abstract 



In the framework of chaotic scattering we analyze passive tracer transport 
in finite systems. In particular, we study models with open streamlines and 
a finite number of recirculation zones. In the non trivial case with a small 



o 

m 

' number of recirculation zones a description by mean of asymptotic quanti- 

^ ' ties (such as the eddy diffusivity) is not appropriate. The non asymptotic 

■ properties of dispersion are characterized by means of the exit time statistics, 
which shows strong sensitivity on initial conditions. This yields a probability 

\ distribution function with long tails, making impossible a characterization in 

O ' terms of a unique typical exit time. 

^ I. INTRODUCTION 

^ \ The problem of transport in velocity fields characterized by different flow regimes in 

c3 \ different subareas (i.e. a "spatially disordered set of streamlines" [^) has been considered by 
many authors. Particular attention has been devoted to steady and time-dependent oceanic 
and atmospheric flows with recirculations and the related problem of the dispersion in 
porous media In this paper we will focus our attention on flows with recirculations 

of geophysical interest. The study of transport properties in presence of recirculations has 
a crucial relevance, since gyre- or eddy-like recirculating patterns are ubiquitous features in 
different areas of the world ocean and atmosphere. In the ocean these features are typically 
induced by forcing spatial structures at the boundaries (e.g. bottom topography, wind stress 
curl, coast-lines) or by intrinsic dynamical reasons (mesoscale eddies), one interacting with 
the other (for a general reference see [^]). 

It is worth noting that large-scale meandering jets, which are typically associated with 
the extensions of western boundary currents, often separate ocean regions characterized 
by different physical and biogeochemical characteristics. Consequently, the study of mixing 
processes in correspondence of them is important also for multidisciplinary investigations |T^ 
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(e.g. the biological effects of longitudinal transport in western boundary current extensions 
see 
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The systems we consider are characterized by the joint presence of open streamline areas 
where particle motion is essentially a ballistic flight and closed streamline regions typically 
distributed according to a periodic geometry, where particles tend to be trapped. The easiest 
way to study the transport properties in such systems is by averaging over smaller space 
or over shorter time scales. Typically, this results in the possibility of describing dispersion 
by an equation for macroscopic quantities of the system such as the average passive scalar 
concentration in terms of an effective drift and an effective diffusivity (as classically done 
by Taylor |jl2[; see |jl]], for a detailed discussion), at least in the case of standard diffusion. 
However, a description in terms of equations for macroscopic quantities needs to average 
out the small scale (fast) features of the velocity field and thus such approach applies only 
for asymptotic times, when particles have been able to thoroughly sample the different flow 
regimes in the system. 

As stressed by Young [|l|, before reaching this asymptotic regime very interesting tran- 
sient behaviors, which cannot be described within an effective diffusion model, could occur. 
The transient regime may be very long and, if the system is finite, the asymptotic one may 
not be reached. This is, typically, the case of finite domain systems with no large scale 
separation between the size of the domain and the largest characteristic Eulerian length. In 
realistic flows, which are usually characterized by a fairly limited number of recirculations, 
fluid particles ordinarily sample just a fraction of the available regimes. 

Therefore, it is often not possible to characterize dispersion simply in terms of asymptotic 
quantities such as average velocity and diffusion coefficients: different approaches are needed, 
as done in [13-1^ for the characterization of transport in closed domains; such as the 



symbolic dynamics approach to the sub-diffusive behavior in a stochastic layer and to mixing 
in meandering jets respectively described in [|16| and |T^; or the study of tracer dynamics 
in open flows in terms of chaotic scattering by Tel and coworkers [^^-|2D|; and the exit time 
approach of 0. 

The aim of this paper is to describe dispersion in finite size systems. In particular, we 
want to characterize flows with a small number of recirculations using some ideas stemming 
from the chaotic scattering theory. 

In Section |I| we briefly review some methods for the study of transport in non asymptotic 
regimes, namely the finite-size diffusion coefficients and the exit time statistics originated 
from the chaotic scattering phenomenon. Section ^TT] contains a description of the two studied 
models (traveling wave and meandering jet) and some numerical results. Section |^ is 
devoted to the comparison between numerical results and a probabilistic model. Conclusions 
are presented in Section 0. 



II. TOOLS FOR THE STUDY OF NON- ASYMPTOTIC 
TRANSPORT PROPERTIES 

The investigation of passive tracer diffusion is usually reduced to the study of an effec- 
tive equation describing the long-time, large-distance average tracer concentration behavior. 
Under rather general hypotheses, given an Eulerian velocity field the long-time transport 
process is uniquely characterized by the effective diffusion (or diffusivity) tensor D^ij-. 
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D^ij = lim — < {xi{t)- < Xi >){xj{t)~ < Xj >) > , (1) 

where x(t) is the position of the tracer particle at time t\ i,j = 1, ■■■,(!, and d is the 
spatial dimension; the average is taken over the tracer initial conditions or, equivalently, 
over an ensemble of tracer particles. The effective diffusion tensor Dfj takes into account 
the molecular diffusivity and the details of the velocity field. Even in presence of simple 
Eulerian fields (e.g. laminar and periodic in time) the diffusion coefficient as a function of 
the parameters of the velocity field can display a rather non trivial behavior PT|J2^ 



It is worth underlining that the diffusivity tensor (Q) is mathematically well defined only 
in the asymptotic limit, therefore its use in finite size domains yields meaningful results 
only if the characteristic length of the velocity field is much smaller than the size of the 
domain. If this is not the case (e.g. in many geophysical settings or plasma physics [|T3[] ), 
dispersion can be characterized more satisfactorily using concepts and techniques borrowed 
by the dynamical systems theory which will be discussed in the following. 

As an example, consider the relative diffusion of a cloud of A/" test particles in a smooth 
velocity field with characteristic length and assume that the Lagrangian motion is chaotic. 
Denoting with R{t) the radius of the cloud: 

R'{t) = ^ j:\Mt) -{At))\\ (2) 
•'^ k=i 



with 

M^)) = ^E^n(t) (3) 



1 ^ 



we expect the following regimes to hold 



ri/2 



^ . i?2(o) exp(L(2)t) if i?2(t) « 



where L{2) > 2A is the generalized Lyapunov exponent is the diffusion coefficient 

under the assumption of isotropy, D^ij = SijD^, and the over-bar denotes the average over 
initial conditions of the cloud. Equation (Q) states that as long as the size of the cloud is 
relatively small, chaotic behavior prevails; later on, dispersion displays a diffusive behavior. 

Another way to describe the above behavior, even though at first glance rather artificial, 
is by introducing the "doubling time" T[6) at scale 6 as follows: define a series of thresholds 
^{n) _ j,ri^{o)^ where S^^^ is the initial size of the cloud, and then measure the time T((5(°)) 
it takes for the growth from 6^^^ to 6^^^ = r5^^\ and so on for T(6^^^) , T{6^'^^) , ... up to 
the largest considered scale. Though strictly speaking the term "doubling time" refers to 
the threshold rate r = 2, any value can be chosen for r, even if a too large one might not 
separate different scale contributions. 

Performing Ai ^ 1 experiments with different initial conditions for the cloud, we define 
an average doubling time T(5) at scale 6 as 

1 M 

T{5) =< T(5) >M= T7 E Tm{S) . (5) 

m=l 
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It is worth noting that the average in (^) is different from the usual time average (see ||24 
for a detailed discussion of this point). A finite size diffusion coefficient D{6) can then be 
introduced as 

D{S) = S'XiS) , (6) 

where 

In r 

MS) ^ ^ (7) 



is the finite size Lagrangian Lyapunov exponent [|T^|. It can be shown that the Lyapunov 
exponent A can be obtained from X{6) for 6^0, namely X{6) = A for 5 -C /«. 

For a tracer cloud of non-infinitesimal size T{6) depends on the details of the nonlinear 
mechanisms of expansion: in the case of standard diffusion D{6) is a constant, i.e. 1/T(5) ~ 
0. Thus 



The fixed scale analysis allows us to extract physical information at different spatial 
scales avoiding some unpleasant consequences resulting from working at a fixed delay time 
t. For instance, in presence of strong intermittency, R'^{t) as a function of t can be rather 
different from one realization to another. In figure la we present an example taken from 



1 14 1 where different exponential rates of growth for different realizations of R (t) produce a 



spurious behavior of R'^{t). Since at small scales, say one typically has an exponential 
growth of R{t) due to Lagrangian chaos and at large scale a diffusive behavior, one can 
mimic R{t) as follows: 

R2,,^ _ / Sy^' R{t) < In 

^ ~ \ 2D{t - Q R{t) > lu ^ ^ 

where is such that R{t) is continuous. 

Figure lb shows the average R^{t) versus time t; at large times the diffusive behavior is 
recovered, but at intermediate times an apparently anomalous regime occurs. This last 
regime is only due to the superposition of exponential and diffusive contributions by different 
samples at the same time. On the other hand the doubling time analysis yields unambiguous 
results, see figure Ic. 

Let us remark that the above technique recovers the usual asymptotic description when 
there is a large scales separation and also if a genuine anomalous diffusion occurs |25], in 
addition, it constitutes a systematic method to treat situations in which the scales are not 
well separated. Moreover, the finite size diffusion coefficient D{6) enables the understanding 
of the different spreading mechanisms at different scales. This has been recently shown 
in [jl5[, where this method is applied to analyzing experimental trajectories described by 
surface drifters in the Adriatic Sea. 

Another interesting approach to the study of transport properties in finite systems is 
the chaotic scattering theory used in |lT8|- pO|] for passive tracer advection in open flows. In 
a nutshell, chaotic scattering can be summarized as follows. A particle arriving from, say. 
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X = —oo enters a region (defined as the scattering, or interacting, region) where due to the 
presence of a potential it scatters, then exits and goes to x = oo. Typically, for a rather 
general class of potentials ||26| the time a particle takes to escape from the scattering region 



can be very sensitive on the impact parameter b and thus displaying a chaotic character 
(like a ball in a pinball game). This justifies the definition in terms of chaos, even if since 
chaos is a time- asymptotic concept, from a technical point of view this kind of behavior is 
not chaotic: particles after a transient (even if very long) leave the scattering region and 
enter a regular motion regime. 

A nice example of the application of chaotic scattering theory in passive tracer study is 
given in [0,11^] for the motion of Lagrangian tracers in blinking vortex-sink system and in a 
von Karman vortex street behind a cylinder in a channel. Tracer particles can temporarily 
be trapped in certain regions, e.g. the wakes of the von Karman street, performing very 
irregular paths. On the other hand, since the non stationarity of the flow is mainly restricted 
to a finite mixing region around the obstacles, the asymptotic almost free particle motion is 
also recovered. 

For the use of the exit time approach for the transport and mixing in volume preserving 



maps see Ref. p7[| . 



The analogy between chaotic scattering, occurring in Hamiltonian systems, and passive 
scalar motion can be drawn in formal terms for two dimensional incompressible velocity 
field |]18|. In this case the Eulerian field is described by a stream-function ip^x, y, t), and the 



corresponding equations for the Lagrangian evolution are: 

dx dip{x,y,t) dy dip{x,y,t) 
dt dy ' dt dx 



(10) 



The equations ([10|) are nothing but the canonical equations for a one dimensional time 
dependent Hamiltonian system, where the stream function plays the role of the Hamiltonian. 

Since in chaotic scattering the particle exit time from the interacting region depends 
strongly on the initial conditions, it is interesting to look at the time delay function [ |18[| , 
i.e. the exit time as a function of the initial position, r(x(0), ?/(0)), e.g. with x{0) = xq. 
In presence of chaotic scattering r(xo,?/(0)), displays a rather irregular shape (see, e.g., fig. 
5.18 in O). 

As already remarked, since in chaotic scattering the irregular character is confined both 
in space and time (i.e. one has the so-called transient chaos), the Lyapunov exponent is 
trivially zero; however, a sort of high sensitivity on initial conditions is suggested by the 
occurrence of very different time delays for very close deployment locations. 

The presence of large excursions for the exit time has an obvious relevance for transport 
processes in finite size systems. The wild variations of T{xo,y{0)) pose severe limits on 
the possibility to make prediction on the particles behavior and force us to use statistical 
approaches. This leads to introduce the probability distribution function, P(r), of the exit 
times t{xo, y{0)). 



III. TWO SIMPLE FLOWS 

In this paper two idealized flow models are studied; both are reminiscent of oceanographic 
features. We consider finite amplitude Rossby waves in a channel, recently revisited in 
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(see also [p9| , p0| ). In the reference frame moving with the phase speed of the wave, 
the flow field shows a central open streamline region (ballistic motion) sided by trapping 
recirculations (see fig. 2a); this flow pattern is a suitable, even if simplified, prototype for 
studying the effect of a finite number of trapping areas on the longitudinal dispersion of 
particles. 



The two-dimensional incompressible Rossby wave flow p8| here considered is specified 
by the following stream-function: 



"^oix, y) = Aq sm{Kox) sin(Loy) - Coy 



(11) 



where Aq is related to the maximum velocity in the direction, {Kq, Lq) is the wave vector 
and Co is the phase speed of the primary wave in the x— direction. In ( pH]) \E'o is expressed 
in the reference frame co-moving with the primary wave. 

In order to reproduce the instabilities usually present in geophysical flows, we introduce 
a time-periodic perturbation, 6'^{x,y,t): 



5\E'(x, y,t) = a sin(_ft'ix — Qt) sin(Liy) 



(12) 



where a > is a (not necessarily small) parameter which controls the amplitude of the 
perturbation and {Ki,Li) is the wave- vector of the perturbation (secondary wave). Even 
though realistic disturbances cannot be characterized in terms of a single wave alone [^], the 
perturbation (|12D is a first step towards a description of the complex structure of transport 
mechanisms in these systems. 

The second flow we investigate is a meandering jet, which represents a natural extension 
of the above Rossby wave system, and it was extensively studied in the hterature with 
particular reference to the Gulf Stream |T0| , |32| . Again, the flow can be subdivided into 
different regions roughly corresponding to a prograde flow (in reality, the current jet core; 
in our schematization, the open streamline regime), recirculation regions and, at a farther 
distance, an essentially quiescent (far) field (see fig. 2b). 

We consider now fluid particle trajectories in the two-dimensional kinematic model origi- 
nally proposed by Bower and thereafter widely studied [p!7|,^^-|36| . The large-scale flow, 
in a reference frame moving eastward with a velocity coinciding with the meander phase 
speed and suitably nondimensionalized, is expressed by the stream function: 



'4^ix,y) 



tanh 



y — B cos kx 
;i + A;2E2sin2 fca;)i/2 



+ cy. 



(13) 



In fig. 2b we show the streamlines in the co-moving frame. 

As mentioned above, chaotic advection may be introduced trough a time dependence. 



Among the different mechanisms discussed in , we chose here a time-periodic oscillation 
of the meander amplitude: 



B{t) = 5o + 7 ■ cos{ujt + 



(14) 



where we set Bq = 1.2, 7 = 0.3, uj = 0.4 and 6 = 7r/2. 
motivated by observations and numerical results, as discussed in 

It is worth underlining that, even if quite a great deal of efforts has been devoted to 



The parameters choice is mainly 



study fluid exchange across the jet [|33|-p^, very little is known as to tracer behavior in the 



6 



along-jet direction (periodic flows with open streamlines have been proposed as model for 
the meandering jets; the presence of recirculations has been seen to induce non-trivial effects 
on the along-jet dispersion, see p7[] ). 

Tracer particle trajectories have been numerically generated from eqs. (0) with the 
stream-functions corresponding to the traveling wave (|lT]-^) and the meandering jet 



14). However, since the results relative to the two flows are qualitatively the same, we shall 



present and discuss just those obtained for the meandering jet. 

As can be seen from the streamfunctions, the two flows are periodic in the longitudinal 
direction. Since we are interested in the characterization of longitudinal transport, the 
number of elementary flow structures (or cells), Nc, constituting the system, is a crucial 
parameter. For very large Nc, the dispersion properties of the system can be studied using 
asymptotic techniques, e.g. the multiscale method |3^J3^. On the contrary, we mainly 
concentrate on systems with a small number of cells, namely — 2 — 10. 

The first focus of our analysis is the particle exit time (or time delay function, see Sect. ^ 
as a function of the initial position, i.e. the time r(xo,?/(0)) a tracer particle deployed at 
{xo,y{0)) takes to reach the boundary Xmax = Nc2TT/k. Two very different scenarios occur 
for large and small Nc. 

Figs. 3a-d show the behavior of r(xo,?/(0)) for Nc = 3, 10, 100, 1000. Increasing the sys- 
tem size there is a clear change in the shape of t{xo, 2/(0)): highly inhomogeneous structures 
(fractal objects) for low Nc (figs. 3a,b), with very strong fluctuations of the exit time value 
even for small variations of the initial conditions. As Nc increases, the shape of r(xo,2/(0)) 
becomes more and more homogeneous. 

The fractal character of t(xo, 2/(0)) is evident from figs. 3a and the enlargements figs. 4a-b, 
which suggest the self-similarity of the structures at different scales. This can be quantita- 
tively assessed studying the correlation dimension V of the initial condition set {y{0)} such 
that T{xo,y{0)) > O. Using the Grassberger and Procaccia algorithm |^0|, i.e. computing 
the percentage C(r) of pairs (y,, yj) such that | yi — yj \< r; for small r we obtain C(r) ~ 
(shown in fig. 5) with T> < 1. The value of T> can depend weakly on the threshold 6, e.g. 
for G = 15 and the parameters of fig. 3a T> results 0.83 and 6 = 30 yields T> = 0.78 

Let us remark that also for very large values of Nc some fractal structures may be 
present, in particular, on very small scales (detectable only for infinitesimally close particles). 
However, we do not consider this feature because on those scales in real fluids we expect the 
presence of smoothing due to molecular diffusion. 

In order to characterize the system behavior, we have also studied the probability density 
function Pn^{t). In figs. 6a-d the probability density functions Pnc{t) corresponding to 
Nc = 3, 10, 100, 1000 are shown. For large Nc (fig. 6c, d) -P/v^(t) displays an asymptotic shape 
which can be obtained with simple probabilistic arguments (see below. Sect. fV\) , whereas 
in the opposite case (fig. 6a,b) the distributions exhibit sharp peaks in correspondence of 
the ballistic time and exponential tails indicating the possibility of very large excursions. 

The above results show that the dispersion process in a finite size system cannot be 
described in terms of a unique characteristic time. As shown in figs. 3a,b, for A^^c = 3 — 10 
the exit time r(xo,?/(0)) exhibits very strong fluctuations within 3 — 4 orders of magnitude, 
which prevents the possibility to extract meaningful information just from the average of 
T{xo,y{0)). This is reflected in the exit times probability distribution: for large Nc the 
shape of the distribution suggests the appropriateness of the average exit time to describe 
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the diffusion, whereas this is not the case for small Nc, where PhX"^) exhibits exponential 
tails. 



IV. A PROBABILISTIC MODEL 
COMPARED WITH NUMERICAL RESULTS 



In the previous Section we have discussed some statistical properties of the Lagrangian 



dynamics generated by the streamfunctions ([TT] , p!3D ; it seems thus natural to look for prob- 
abilistic models reproducing the above properties. 

The Rossby wave flow (fig. 2a) shows two different regions; the central one characterized 
by ballistic motion (open streamlines) and particle trapping recirculations (closed stream- 
lines) on both sides of the ballistic regime. The meandering jet flow (fig. 2b), in addition to 
the recirculation and ballistic ones, presents two far field regions moving retrogradely with 
respect to the jet core. However, with our parameter choice the far field is practically never 
visited by the tracer particles Jl^^ill and thus we shall not consider them. 



This intrinsic subdivision of the flow fields suggests to build a discrete symbolic picture 
of the motion: one can label the trapping and ballistic regions with the symbols T and 
B respectively and then construct a probabilistic model. Since the studied flows are time 
periodic we express the time in number periods of the perturbation, namely ^-n jio. 

The simplest conceivable probabilistic model is a Bernoulli scheme such that at each time 
the particle can be in the states B with probability p or T with probability 1 — p. However, 
we do not expect such a model to give a good description of the system, because once the 
length L of the system (i.e. the equivalent of is fixed, the probability PlIt) that at time 
T > L the particle exits from the domain depends just on the parameter p, and just one free 
parameter is obviously not enough to characterize the system behavior. 

The natural choice for a model which take into account memory effects is a Markov 
chain [^. The process is completely defined by the transition matrix Wij, i.e. by the 
probability to go in one step to the state j starting from the state i = B, T) which has 
the properties: 

Wij > and J2 = 1 • (15) 

j 

The probabilities Pj to be in state i are given by 

P, = Y.P,W,,. (16) 

j 

Now we consider the following stochastic process for the evolution of a particle 

^ ^ [0 with probability WiT , 

where i represents the state visited at time t — 1 and Ax{t) the increment in the position of 
the particle at time t. 

Then we can compute the probability Pbij) as follows: 
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r-1 



PLir)= E PL-iik) FBBir - k) (18) 

k=L-l 

where FBsin) is the probabihty of first arrival from state B to the state B inn steps, which 



obeys the following recursive formula |^ 



n-l 



FBB{n) = {W^)bb - E PBsin - k) iW')BB , (19) 



k=l 



where W'' indicates the A;-th power of the matrix W. Applying recursively the ( [T^ ) yields 
PLir)= E PBsir-k,) e' ^bb(A:i-A;2)... E PBBikL-2 ~ k^.,) P,{kL-i) (20) 

fci=L-l fe2=-L-2 fcL|=l 



being 



^^^^^ ~ { Pt{Wtt)''-'Wtb k>2. ^^^^ 



Even if, in general a one-step Markov process is not enough for a detailed description 
of statistical properties of the dynamical system |T7| , p2|^ , it may result a good description 



as long as r is not too small: one indeed expects that for large exit times the memory 
effects are less important. The expression (|20|) for -Pl(t) has been compared with the exit 
time probabilities numerically computed in Sect. |IT|. In order to carry out the comparison 
between the probabilistic model and the numerical results we have first to evaluate the 
parameters of the model, i.e. L, Pi, Wij. 

In the numerical evaluation of the matrix Wij and the probabilities Pi we have proceeded 
according to the following scheme. Expressing the time in number of periods of the per- 
turbation, the transition probabilities are computed from a long trajectory xo,xi,---,x„ 
(n ^ 1) in an infinite system {Nc = oo) as 

W,, = hm (22) 

where Nn{i) is the number of times that, along the trajectory, the particle visits the state 
i {i = T OT B) and Nn{i,j) is the number of times that x^ is in state i and Xf+i is in 
state j = T,B). The visit probabilities Pi are simply given by lim„^oo -^n(0/''^- The 
identification of the visited state is performed by controlling the value of the stream function 
and the sign of the velocity along the x, y-axis fl^ . 

In order to evaluate L we have to take into account that in the physical system the 
spatial increment per time step (i.e. during one period of the perturbation). Ax may vary 
with time. We have then computed the most probable value Ax, of Xt+(2n/uj) — Xt, hence we 
have rescaled L using [L/Ax] = L, the square brackets [■] here indicate the integer part of 
the argument. 

In figs. 6a-b the probability density distributions Pl{t) are compared with those obtained 
using (BO). The small shift of the peak at the beginning of the distributions could depend 
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on the estimate of Ax. As can be seen from the figures 6a-b, even for L ~ 3 — 10 we obtain 
a good description of Piij) at least for r ^ 1, i.e. for those particles that experience a 
large number of transitions between different states. Increasing L, Pl{t) is better and better 
approximated by the Markov chain prediction. It is worth remarking that we are comparing 
the numerical results with the prediction of the probabilistic model without performing any 
data fitting. 

On space/time scales much larger than the typical time/ space scales of the velocity field 
the evolution of a test particle follows a diffusive scenario. Thus we expect the following 
process to be a good model for our systems at least in the limit L ^ 1 and r ^ 1: consider 
the diffusive stochastic process 

x{t) -x{0) =vt + V2Dw{t) (23) 

where 

_ ^ l^ {x{t) -x{0)) ^ ^ ^ 1_ ^ ^^^^^ _ ^^^^ _ _^^2 ^ ^24) 

w(t) is a Wiener process, i.e. a Gaussian process with w{0) = and 

{w{t)) = , {w{t)w{t')) = mm[t,t']. (25) 
Defining r the first exit time i.e. the maximum time for which 

x{t) < Xmax =VT + V2Dw{t) (26) 

the probability density of r can be calculated M] obtaining 



p / \ I •^max I ^ 
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[vT Xjjiax ) 



(27) 



The maximum of Px^axi'^) is reached for r^ax which can be estimated, if D/v is not too 
large, as 

_ Xmax /2S^ 

'max — _ • \^'~'J 
V 

The quantities v and D can be evaluated in terms of the exit time statistics as follows 

D = ^ (29) 



< r >e 2 < r > 

where the ( ■ )g indicates the average over the ensemble of Np particles, i.e. = 

In figs. 6a-d the probability distributions calculated in Section |T| are compared, at 
different values of L = Xmax, with the results given by (|27|) with v and D obtained by (|29D 
for large Nc. 

The larger L, the better is the fitting with numerical data, nevertheless the probability 
density function p7| ) does not capture the tail behavior of the physical system probability 
PL(r). This is not very surprising because the process described by eq. (|23|) is 5-correlated 
in time and thus it is not able to describe long-range correlated events responsible of 
such a tail. Moreover, at variance with the physical system, in model ( p3D the velocity does 
not have any bound. 
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V. CONCLUSION AND DISCUSSIONS 



In this paper we have studied non asymptotic properties for passive tracer transport. For 
a finite system (of size L), which is a rather common case in real problems (e.g. geophysical 
and plasma fiows), the usual characterization by means of asymptotic quantities (such as 
the eddy diffusivity) is not appropriate if L is not very large with respect to the typical 
length scale of the velocity field. 

We have considered two models of geophysical interest (traveling waves and meander- 
ing jet) studying their Lagrangian transport properties at varying the longitudinal size L. 
Transport has been analyzed in terms of the statistical properties of the exit times of parti- 
cles from the systems. This analysis has been carried out borrowing concepts from chaotic 
scattering theory. 

In the limit of very large L the usual asymptotic scenario is recovered, i.e. the mean 
velocity v and the diffusion coefficient D completely characterize the transport process. In 
this case, one typical time is enough to describe the basic features of the process. On 
the contrary, in the more interesting (and realistic) cases with a not very large L, there 
is no unique characteristic time. Indeed the exit times show a strong sensitivity to initial 
conditions (this is a manifestation of transient chaos in a non-chaotic system), limiting the 
possibility of a detailed forecasting of particle behavior. This limitation is mirrored in the 
probability distribution function of the exit times which displays rather long tails, making 
impossible the characterization in terms of a unique time (i.e. the average exit time). 

Suitable processes (e.g. Markov chains) prove to capture the relevant statistical aspects of 
transport process. If L is very large the exit times statistics is well described in terms of first 
exit time problem for a linear Langevin equation involving only v and D. For systems with 
a moderate number of recirculation zones one has to introduce a more detailed probabilistic 
model. 
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TABLES 



Table 1 

Transition matrix elements 



Wbb 


.66 


Wbt 


.34 


Wtb 


.12 


Wtt 


.88 



Table 2 

Visit probabilities 



Pb 


.26 


Pt 


.74 



The probabilities are evaluated for the stream function ( |I3D with parameters k = 
An/15, Bo = 1.2, c = 0.12, uj = 0.4 and 7 = 0.3. 

The statistics have been computed over 2 10^ periods. 
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FIGURE CAPTIONS 



FIGURE 1 :(a) Three realizations of R^{t) built as eq. (|) with: 7 = 0.08,0.05,0.3 and 
5q = KrJ, D = 1.5. 

(b) -R^(t) as function of t averaged on the three realizations shown in figure la. The 
apparent anomalous regime and the diffusive one are shown. 

(c) X{6) vs 6, with Lyapunov and diffusive regimes. 

FIGURE 2 :(a) Streamlines for the time- independent Rossby wave flow with: Aq = Lq = 
Kq = 1, Co = 0.5, a = 0. The labels T and B refer to the trapping regions (recircula- 
tions) and ballistic ones (the jet channel) 

(b) The same as figure (a) for the meandering jet flow with: k = 47r/15, Bq = 1.2, c = 
0.12 and 7 = 0. 

FIGURE 3 r(xo,?/(0)) for the meandering jet flow with parameters for the unperturbed 
streamfunction as fig.2b and u = 0.4, 7 = 0.3 9 = n/2, and xq = 0.1, y{0) G [—3,3] 
for (a) Ne = 3 cells, (b) = 10, (c) A^, = 100 and (d) (a) = 100. 

FIGURE 4 Two enlargements of figure 3a: (a) for y{0) G [1.15 : 1.20] and (b) y{0) G 
[1.1910 : 1.1925]. 

FIGURE 5 Correlation integral C{r) versus r computed from the data of fig. 3a, with 
threshold B = 15. The straight line has slope V = 0.83. 

FIGURE 6 Pn^{t) versus r for the meandering jet flow with the same parameters as in 
fig. 3, computed with 410"^ particles starting in x(0) G [0,7r/fc], y{0) G [—3,3] for (a) 
Nc = 3 cells, (b) A^c = 10, (c) = 100 and (d) = 1000. The prediction (^) is 
shown as dashed continuous line in all the cases. The parameters v and D have been 
calculated from (p9|) for Nc = 1000 and evaluated as 0.365 and 1.038 respectively. In 
figure (a) and (b) the Markovian prediction is also shown as continuous lines using 
Ax = 1.36 and evaluating the probabilities Pi and Wij as reported in tables 1 and 2. 
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